Gaussian-state quantum-illumination receivers for target detection 
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The signal half of an entangled twin-beam, generated using spontaneous parametric downcon- 
version, interrogates a region of space that is suspected of containing a target, and has high loss 
and high (entanglement-breaking) background noise. A joint measurement is performed on the re- 
turned light and the idler beam that was retained at the transmitter. An optimal quantum receiver, 
whose implementation is not yet known, was shown to achieve 6dB gain in the error-probability 
exponent relative to that achieved with a single coherent-state (classical) laser transmitter and the 
optimum receiver. We present two structured optical receivers that achieve up to 3dB gain in 
the error exponent over that attained with the classical sensor. These are to our knowledge the 
first designs of quantum-optical sensors for target detection, which can be readily implemented in a 
proof-of-concept experiment, that appreciably outperform the best classical sensor in the low-signal- 
brightness, high-loss and high-noise operating regime. 

PACS numbers: 



A distant region engulfed in bright thermal light, sus- 
pected of containing a weakly reflecting target, is inter- 
rogated using an optical transmitter. The return light 
is processed by a receiver to decide whether or not the 
target is present. Recent work [TJ|2J[3] has shown that in 
the above scenario, a "quantum illumination" (QI) trans- 
mitter, i.e., one that generates entangled Gaussian-state 
light via continuous-wave pumped spontaneous paramet- 
ric downconversion (SPDC), in conjunction with the 
optimal quantum receiver, substantially outperforms a 
coherent-state (un-entangled) transmitter and the cor- 
responding optimum-measurement receiver. This advan- 
tage accrues despite the loss of entanglement between the 
target-return and the idler beams due to the high loss and 
noise in the intervening medium. This is the first example 
of an entanglement-based performance gain in a bosonic 
channel where the initial entanglement does not survive 
the loss and noise in the system. The SPDC transmit- 
ter and optimal receiver combination has been shown to 
yield up to a factor of 4 (i.e., 6 dB) gain in the error- 
probability exponent over a coherent state transmitter 
and optimal receiver combination, in a highly lossy and 
noisy scenario [3 . The optimal receiver for the former 
source corresponds to the Helstrom minimum probabil- 
ity of error (MPE) measurement [3] under two hypotheses 
- Ho: target absent, and H\: target present. It can be 
expressed as a projective measurement onto the positive 
eigenspace of the difference of the joint target-return and 
idler density operators under the two hypotheses. How- 
ever, no known structured optical receiver is yet able to 
attain the full 6 dB predicted performance gain. 

In this paper we present two structured receivers, 
which, when used in conjunction with the SPDC trans- 
mitter, are shown to achieve up to a factor of 2 error- 
exponent advantage i.e., half of the full factor of 4 
predicted by the Helstrom bound — over the optimum- 
reception classical sensor, in the low signal brightness, 
high loss and high noise regime. The first receiver uses a 
low-gain optical parametric amplifier (OPA) and ideal 



photon counting [5], whereas the second uses phase- 
conjugation followed by balanced dual detection. Both 
receivers attain the same asymptotic error exponent, al- 
though the second receiver yields slightly better perfor- 
mance than the first. Both receivers attempt to detect 
the remnant phase-sensitive cross correlation between the 
return-idler mode pairs when the target is present [6]. 
Both of our proposed receivers, consisting of separa- 
ble measurements over M pairs of target-return and 
idler modes, offer strictly better performance than any 
classical-state transceiver, and have low-complexity im- 
plementations. Apart from the binary-hypothesis target- 
detection problem considered here, our receivers have 
been shown to considerably outperform conventional op- 
tical sensing and communications applications built on 
the quantum illumination concept, such as two-way se- 
cure communications [7] and standoff one-vs-two-target 
resolution sensing [8]. 

Consider M independent signal-idler mode pairs ob- 
tained from SPDC, {effi , a^}; 1 < k < M. Each T- 
sec-long transmission comprises M = WT 3> 1 signal- 
idler mode pairs, where W is the SPDC sources phase- 
matching bandwidth. Each mode pair is in an identi- 
cal entangled two-mode-squeezed state with a Fock-basis 
representation [9 
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where N$ is the mean photon number in each signal 
and idler mode. \ip)si is a pure maximally-entangled 
zero-mean Gaussian state with covariance matrix V SI — 
([as ai a s a\) T [a s a\ a s a/ ]) given by 

N s + 1 y/N s (N s + l)~ 

jV s + l y/N s (N s + l) 
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Under hypothesis Hq (no target), the target-return mode 
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&rt = o-Bi where &b is in a thermal state with mean 
photon number Nb ^> 1. Under hypothesis Hi (tar- 
get present), a R — yfnas + \/l — kclb, where k -C 1, 
and as is in a thermal state with mean photon num- 
ber Nb/(1 — k), such that the mean noise photon num- 
ber is equal under both hypotheses. Under Hi, each of 
the M return-idler mode pairs are in a zero-mean Gaus- 
sian state, Pbj, with the covariance matrix of each mode, 
V RI = {[a R a j a R a\] T [a R a\ a R a T ]), given by 



kN s +N b + 1 ^/kN s (N s + 1) 

N s + 1 y/nNsfNs + l) 

y/K.N s (N s + l) kNs+N b 
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Under Hq, the joint return-idler state for each of the M 
mode pairs, p R j, is a product of two zero- mean ther- 
mal states (pn b ® Pn s ) with mean photon numbers Nb 
and respectively, viz., V RI = diag(A^s + l,Ns + 
l,N B ,N s ). 

The binary detection problem is the MPE discrimina- 
tion between Hq and Hi using the optimal measurement 
on the M return-idler mode pairs, (p R j > )^' M for m = 
or 1. The minimum probability of error is given by, 



(M) 



[1— X) n 1/2, where 7„ are the non-negative 
eigenvalues of (/3^) 0M - {Pr])® M % The quantum 
Chernoff bound (QCB), given by Qqcb — mino< s <i Q s 
where Q s = Tr[(/5^]) s (/5^]) 1_s ] , is an upper bound to 

min ano - 1S asymptotically tight in the exponent of the 
minimum error probability [10 . In particular, we have 
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FIG. 1: (Color online) The figure shows five plots of error- 
probabilities and bounds thereof as a function of M. The 
two solid curves marked by arrows (from top to bottom re- 
spectively) are the coherent-state (blue) and the entangled 
Gaussian-state (red) transmitters. The third solid curve 
(black, in between the aforementioned QCB curves) plots the 
error-probability performance of the OPA receiver, whereas 
the dash-dotted curve shows the performance of the phase- 
conjugate receiver. The curve plotted with circles depicts 
the error-probability performance of the coherent-state trans- 
mitter and homodyne detection receiver, which is in fact a 
lower bound to the performance of an arbitrary classical-state 
transmitter, including classically correlated signal-idler trans- 
mitter states. The parameters used to generate the plots are 
N s = 0.01, N B = 20, and k = 0.01. 



where the first inequality (QCB) is asymptotically tight 
as M — > oo [TS]. The QCB is customarily represented 

as -Pg^5; n < e~ MBjQ /2 in terms of an error exponent 

Rq = — ln(QQCB)- The second inequality is a looser up- 
per bound known as the Bhattacharyya bound. Symplec- 
tic decomposition of Gaussian-state covariance matrices 
was used to compute the QCB 0IH], and it was shown 
that in the high loss, weak transmission and bright back- 
ground regime, i.e., with Ns <C 1, K *C 1, and Nb 3> 1, 
the entangled transmitter yields a QCB error-exponent 
Rq = kNs/Nb, which is four times (or 6 dB) higher than 
the error-exponent Rc — kNs/{^Nb) for a coherent- 
state transmitter with a mean photon number Ns per 
mode. In Fig. [T] we plot the QCB for the entangled and 
coherent state transmitters, showing a clear advantage of 
quantum over classical illumination. 

When a coherent-state transmitter is used, each re- 
ceived mode aji is in a thermal state with mean pho- 
ton number Nb, and a mean- field (a R ) = or ^kN s 
for hypotheses H and Hi respectively. Homodyne de- 
tection on each received mode a* R yields a variance- 
(2Nb + l)/4 Gaussian-distributed random variable Xk 



with mean or \/kNs given the hypothesis. The mini- 
mum error probability decision rule is to compare X — 
Xi + . . . + Xm against a threshold: "H " is declared if 
X < (M^/kNs) /2 and "Hi" otherwise. The correspond- 
ing probability of error is 

(Af) 1 f / / kN s M \ e - MRc »™ 
e>om 2 erK W B + 2 J ~ 2^vMR Ch J 

where erfc(:r) = (2/y/n) e~' 2 di, i?c hom = 
kNs/(4:Nb+2) is the error exponent, and the ap- 
proximation holds for kN s M/(4:Nb + 2) > 1. When 
Nb ^S> 1, Rc hom ~ kNs/4:Nb — Rc, so mode- by- mode 
homodyne detection is asymptotically optimal for the 
coherent-state transmitter. 

A. The OPA Receiver 

Unlike the coherent-state transmitter, the entangled 
transmitter results in zero-mean joint return-idler states 
under both hypotheses. The sole distinguishing factor 
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Optical 
Parametric 
Amplifier 

(OPA) 

Small gain, 
G = 1 + E 2 



Total number of clicks registered 
in M received modes 

\ M 

fc=l 



PN\H {n\H ) an d P/v|ffi("|-ffi) using M i.i.d. observa- 
tions is bounded above by the classical Bhattacharyya 
bound [5], 



g( fc ) _ y/Gaf) + y/G - la% k) where with a gmall QPA gain G = 1 + e 2 , e < 1, 



(5) 



the 



FIG. 2: In the OPA receiver, the return modes and idler 
modes are inputs to an optical parametric amplifier (OPA) 
with gain G. The total number of photons, TV, are counted at 
one output port. The receiver decides in favor of hypotheses 
Ho if N is below a threshold Nth, and in favor of Hi otherwise. 



between the two hypotheses that makes quantum illu- 
mination perform superior to the uncntanglcd coherent- 
state transmitter, are the off-diagonal terms of V RI bear- 
ing the remnant phase-sensitive cross correlations of the 
return-idler mode pairs when the target is present. The 
OPA receiver uses an optical parametric amplifier to com- 



bine the incident return and idler modes a R and a 
1 < k < M, producing output mode-pairs: 



and 



5<*> = y/GhP + VG^laf) 



dW=VGaW + VG=Ta\ (k \ 



(3) 
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where G > 1 is the gain of the OPA (see Fig. [2]). Thus, 
under both hypotheses c'*' is in an independent, iden- 
tical, zero-mean thermal state, p c = YlnLo{Nm/ (1 + 
Nm) 1+n )\n)(n\, for m € {0,1}, where the mean photon 
number is given by iV = QN S + (G - 1)(1 + N B ) un- 
der H , and Ni = GN S + (G - 1)(1 + N B + nN s ) + 
2 V /G(G - l)y/KNs(N s + 1) under H v 

The joint state of the M received modes 
{c (fc) ,l < k < M), is the M-fold tensor product 
pf M , and the M-fold product of thermal states is 
diagonal in the M-fold tensor-product of photon-number 
bases. Therefore, the optimum joint quantum measure- 
ment to distinguish between the two hypotheses is to 
count photons on each output mode and decide 
between the two hypotheses based on the total photon 
count N over all M detected modes, using a threshold 
detector. The probability mass function of N under the 
two hypotheses is given by 



PN\H m {n\H m ) = 



M 
n 



1 



N" 



(l + iV m )»+ M ' 



where n = 0, 1,2, . . . and m = or 1. The mean and 
variance of this distribution are MN m and Mer 2 ^ respec- 
tively, where — N m (N m + 1). The minimum error 
probability to distinguish between the two distributions 



error exponent R B is given by 



2 kN s (N s + 1) 



2N S (N S + 1) + 2e 2 (l + 2JV S )(1 + N s + N B ) 
kN s /2N b , 



(6) 



for a choice of e 2 = Ns/y/W^, for N s < 1, k < 1, 
N B 3> 1 [13 . Therefore by construction, for a weak 
transmitter operating in a highly lossy and noisy regime, 
the OPA receiver achieves at least a 3 df3 gain in error ex- 
ponent over the optimum-receiver classical sensor whose 
QCB error exponent R c = kN s /AN b . For N s < 1 and 
e <C 1, both N and Ni •< 1. Hence, a single-photon 
detector (as opposed to a full photon-counting measure- 
ment) suffices to achieve the performance of the receiver 
depicted in Fig. [2] Due to the central limit theorem, for 
M >• 1, P N \ Hm (n\H m ), m £ {0,1} approach Gaussian 
distributions with mean and variance MN m and Ma^ 
respectively. Hence for M > 1. 



e,OPA 



1 



erfc V RqpaM 



e -AffioPA 

2y/nMRoPA ' 



where an error-exponent i?oPA = {N\ — N ) 2 /2(<j + vi) 2 
can be achieved using a threshold detector that decides 
in favor of hypothesis Ho if N < N t h, and Hi other- 
wise, where 7V th = \M{oiNo + a Ni)/(a a + oi)] . Fig. [l] 
shows that -F^OTA is strictly smaller (by 3 dB in error- 
exponent) than -Pg ifLj — the error probability achieved 
by the coherent state transmitter with a homodyne detec- 
tion receiver. One can show using convexity arguments 
that in the high background regime, P^l m is in fact a 
strict lower bound to the error probability achievable by 
an arbitrary classical-state transmitter, which includes 
classically-correlated signal-idler mixed states (i.e., those 
that admit a Glauber P-representation) . 



B. The Phase-Conjugate Receiver 

The phase-conjugate (PC) receiver is another receiver 
whose error-probability achieves the same 3 dB error- 
exponent gain over the optimal classical transceiver in 
the asymptotic operating regime Ns <C 1, K <C 1, N B ^> 
1, and has slightly better performance than the OPA 
receiver (see Fig.[T|. As illustrated in Fig. 13} the receiver 

phase-conjugates all M return modes a/^ , 1 < k < M 
according to 



s (*) 



V2a^ 



(7) 
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FIG. 3: In the PC receiver, the phase-conjugated return 
modes and the idler modes are inputs to a balanced differ- 
ence detector. If the difference in the total number of clicks, 
N, over all M received modes is less than a threshold N t h, 
the receiver decides Ho, and Hi otherwise. 

where ay" 1 are vacuum-state operators needed to preserve 
the commutator. The conjugated return and the retained 
idler are then detected by a dual, balanced difference 
detector: the output modes of the 50-50 beam splitter, 
af = {af + 4 fc) )/V2 and = {af - af ] )/y/2, are 
detected and fed into a unity-gain difference amplifier, 
such that the final measurement is equivalent to 

where N {k) = a^ k) a (k) and N^ = a^ k) a (k) . The final 
decision is based on the sum of the photon counts N over 
all M modes. 

To simplify the subsequent analysis, let us define 

N x ± (NP), N Y = (fiW), N C = (a^a^), and 
Nj = (a\^af^). Under hypothesis Ho, the modes 

and hj are in product thermal states, whereas 
under Hi they are in a zero-mean joint Gaussian state 
with nonzero phase-insensitive cross correlation given by 
(a^ k) a (k) ) = ^JkN s {N s + 1) = C q . Measurement of 

M k \ 1 < k < M, produces a sequence of i.i.d. random 
variables with mean and variance given by N Q = 
and al = N X (N X + 1) + N Y (N Y +_1) under hypothesis 
H , and N x _= 2C q and a\ = N X (N X + 1) + N Y (N Y + 
1) — (A^ c — A/j) 2 /2_under hypothesis Under hypoth- 
esis H we have N x = Ny — (Nc + Nj)/2, whereas 
^Vx = (Nc + Ni)/2 + C q , and N Y = (N c +_ Ni)/2 - C q 
holds for hypothesis Hi . Finally, we have Nc = 1 + Nb 
for H , Nc = 1 + uNs + Nb for Hi, and because the 
idler is unaffected under either hypothesis, Nj = N$- 
For large M, and hypothesis H m , m £ {0, 1}, the distri- 



bution of N — 53fe=i Nk approaches a Gaussian distribu- 
tion with mean and variance given by MN m and Ma^ 
respectively. Therefore the probability of error 

where an error-exponent i?pcR = (iVj — N ) 2 /2(ao + (Ti) 2 
can be achieved using a threshold detector that decides in 
favor of hypotheses Ho when N < N t h and in favor of Hi 
otherwise, where Nth = \M(aiNo + aoNi)/(a Q + Ci)] . 
The corresponding error-exponent is given by 

= kN s (Ns + 1) 

PCR 2N B + 4:N S N B + 6N S + 4k7V| + 3nN s + 2 
« kN s /2N b , (9) 

where N s < 1, k < 1, N B > 1. 

The PC receiver achieves the same 3 dB error-exponent 
gain as the OPA receiver over the optimum-reception 
classical transceiver, though the performance of the for- 
mer is slightly better in absolute terms (see Fig. [TJ . One 
reason for this is that balanced dual-detection cancels 
the common-mode excess noise in a x and ay , which is re- 
flected by the negative term (Nc — Ni) 2 /2 in the variance 
of Nk under hypothesis Hi. On the other hand, the OPA 
receiver operates at very low gain, thus requires much less 
pump power than unity-gain phase-conjugation. 

In summary, we have proposed two receiver structures, 
both viable for low-complexity proof-of-concept experi- 
mental demonstrations using off-the-shelf optical compo- 
nents, which in conjunction with the SPDC entangled- 
state source, could substantially outperform classical 
transceivers for various entangled-state optical sensing 
applications, such as standoff target detection, one-vs- 
two-target resolution sensing [5J and two-way secure com- 
munications [7]. 
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